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Abstract 

Inclusive production of D*^ mesons in two-photon collisions was measured by 
the L3 experiment at LEP. The data were collected at a centre-of-mass energy 
y/s = 189 GeV with an integrated luminosity of 176.4 pb~^. Differential cross sec- 
tions of the process e"'"e~ —>■ e"'"e~D*^X are determined as functions of the trans- 
verse momentum and pseudorapidity of the D*^ mesons in the kinematic region 
1 GeV < p^* < 5 GeV and \ri^*\ < 1.4. The cross section integrated over this phase 
space domain is measured to be 132±22{stat.)±2Q{syst.) pb. The differential cross 
sections are compared with next-to-leading order perturbative QCD calculations. 
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1 Introduction 



The study of charm production in two-photon colhsions provides a means for testing pertur- 
bative QCD and for probing the gluon content of the photon |Tj. Charmed quarks can be 
produced in "direct-photon" processes, in which the interacting photons behave as point-hke 
particles and couple directly to a charmed quark pair. Another class of processes contributing 
to the charm production are the "resolved-photon" processes, where one or both interacting 
photons fluctuate into a flux of partons. In the "single resolved-photon" processes the unre- 
solved photon interacts with a constituent parton from the resolved photon, whereas in the 
"double resolved-photon" processes a hard scattering between the constituent partons of the 
two resolved photons takes place. In the next-to-leading order QCD only the sum of direct and 
resolved-photon processes is unambiguously defined. The experimental measurement of differ- 
ential cross sections for production of open charmed particles allows a detailed investigation of 
the charm production mechanism. 

Charm production in two-photon collisions has been measured at lower centre-of-mass en- 
ergies at PEP, PETRA, TRISTAN and LEP jSHHj, identifying charmed quarks by detecting 
D*^ mesons, soft pions, inclusive leptons and Kg mesons. In a previous measurement by the 
L3 experiment [Oj, events containing charmed quarks were tagged by detecting electrons and 
muons from semileptonic decays of charmed hadrons. In the present study charmed vector 
mesons D*(2010)^ are identified by the small energy released in D* decay, applying the mass 
difference technique ^01 to the decay chains 

I y K-7T+ (1) 

I . K-7r+7r° (2) 

The presence of a low-momentum, "soft" pion, vTg", ensures that the resolution of the mass 
difference M(D°7rg") — M(D°) is superior to the resolution of the reconstructed D° and D*"*" 
masses themselves. The D*+ signal appears as a narrow peak close to the kinematic threshold in 
the mass difference distributions M(K-7r+7r^) - M(K-7r+) and M(K-7r+7r°7r^) - M(K-7r+7r°). 
The combined branching fractions are BR(D*+ D°7r^) ■ BR(D° -> K-7r+) = 0.0263 ± 0.0008 
and BR(D*+ D°7r^) ■ BR(D° K-tt+tt") = 0.0949 ± 0.0064, as given in Reference Hlj. 

2 Selection of Hadronic Two-Photon Events 

The data were collected by the L3 detector (I2j at LEP in 1998 at a centre-of-mass energy 
■^/s = 189 GeV. The integrated luminosity is 176.4 pb~^. 

For efficiency studies, samples of e"'"e~ — >• e"'"e~7*7* e"'"e~ccX events are generated using 
the PYTHIA ^3] and the JAMVG Monte Carlo generators. The background sources are 
simulated by JAMVG (e+e" ^ e+e-r+r"), KORALZ [H] (e+e" ^ ^+^"(7)), KORALW ^ 
(e+e- ^ W+W" ^ f f'f f) and PYTHIA (e+e" ^ qq(7), e+e" ^ e+e"qq ). The Monte Carlo 
events are processed in the same way as the data. 

Reconstruction of the decay chains (1) and (2) requires a sample of events containing 
hadronic final states. Events of the type e"'"e~ e"'"e~7*7* — e~^e~ hadrons are selected by 

The charge conjugate reactions are included throughout the paper. 
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cuts on the energy measured in the electromagnetic and hadron calorimeters and using track- 
ing information. To exclude annihilation events, the total visible energy must not exceed 0.4 ^/s, 
the energy deposited in the electromagnetic calorimeter must be less than 30 GeV and the en- 
ergy in the hadron calorimeter less than 40 GeV. The transverse component of the missing 
momentum vector must be less than 10 GeV and the value of the event thrust must be smaller 
than 0.95. Events are required to have at least three charged particles reconstructed in the 
tracking chamber. 

A total of 1253890 events pass the hadron selection cuts. The contamination from anni- 
hilation processes and two-photon production of tau pairs is less than 0.5%. The subsequent 
reconstruction, which forms D*+ candidates from three-prong decays with invariant mass ex- 
ceeding 2 GeV, suppresses these background contributions to a negligible level. 

The trigger efficiency for detecting two-photon hadronic final states is (87±3)%, determined 
from the data sample itself using a set of independent triggers. 

3 Mass Reconstruction of D*+ Decays 

The identification of D*+ mesons proceeds through two steps: selection of D*' candidates, which 
are then combined with another track to form D*"*" candidates. 

Tracks are used for reconstruction of D° decays if they satisfy the following requirements: 

• Transverse momentum greater than 150 MeV. 

• At least 40 wire hits measured by the tracking chamber. 

• Distance of closest approach to the event vertex smaller than 1 mm in the transverse 
plane. 

A pair of tracks of opposite charge is required to pass the following criteria in order to be 
considered K 7r+ system from a D° decay: 

• The intersection point of the tracks in the transverse plane must be displaced by no more 
than 3 mm away from the event vertex. 

• Pk ■ -Ptt > 2 ■ 10"'^, where Pk and P^r are the probabilities, calculated from the measured 
energy loss dE/dX of each track, for kaon and pion mass hypotheses of the corresponding 
tracks. 

The selection of tracks and their combinations into neutral pairs is identical for the channels 
(1) and (2) in order to minimize the relative systematic error between the two decay modes. 

To reconstruct D'' decays in the K~7r+7r° decay mode, a neutral pion is added to the selected 
K~7r^ system. Neutral pion candidates are formed by a pair of photons, identified as isolated 
clusters in the electromagnetic calorimeter, not matched with a charged track. Photons are 
accepted for 7r° reconstruction if they are detected in the barrel part of the electromagnetic 
calorimeter and their energies are greater than 100 MeV. The 7r° candidates must have the 
invariant mass of the photon pair in the mass window of ±15 MeV around the vr'' mass. The 
decay D° — »• K~7r+7r° proceeds dominantly through one of the quasi-two-body intermediate 
states K*°7r°, K*^7r+ and K~p^ ^l]. We require either the invariant mass of a Kvr subsystem 
to be within ±80 MeV of the corresponding K*(892) mass or the invariant mass of the 7r"'"7r° 
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system to be within ±150 MeV of the mass. If this condition is met for a given intermediate 
resonant state, we make use of the P-wave properties of a vector particle decay into two scalar 
particles and demand in addition the helicity angle 9* of the corresponding decay cascade to 
satisfy the condition |cos 6**1 > 0.4. The helicity angle 9* is defined as the angle between the 
direction of a decay product of the vector resonance (K**^, K*~ or p+) and the direction of 
the pseudoscalar particle (vr^, vr^ or K~) from the D° decay, calculated in the intermediate 
resonance rest frame. 

To reduce the combinatorial background when reconstructing D° decays into the K~7r"'" 
final state, the opening angle of the track pair in space must be smaller than 2.5 rad. The 
combinatorial background for the K~7r"'"7r° decay mode is suppressed by requiring the solid 
angle, defined by the directions of flight of the three decay products, to be smaller than 2 srad. 

The invariant mass of the K~7r+ system is calculated and if it is in the range of ±100 MeV 
around the mass of the D° mesons [TT], the combination is retained as a D° candidate for the 
decay channel (1). The corresponding mass window for candidates in channel (2) is ±50 MeV. 
The different mass windows reflect the corresponding D'^ mass resolutions, as obtained by Monte 
Carlo studies. The better resolution of the D° reconstruction in channel (2) is due to the softer 
and thus better measured charged particles produced in the three body decay and to the use 
of a well-measured 7r°. 

Finally, the probability that a particular K^7r+ combination comes from a D° decay in 
channel (1) is determined from a IC kinematic fit, in which the invariant mass of the pair is 
constrained to the D° mass. For the K~7r"'"7r° final state we perform a 2C fit, constraining 
the mass of the whole system to the D° mass and the two-photon mass to the 7r° mass. A 
combination is accepted as a D° candidate if the confidence level of the fit is greater than 0.5%. 

In the second step of the D*"*" reconstruction, we consider all combinations of a given D° 
candidate with an additional track of positive charge, assumed to be the soft pion TTg", resulting 
from the D*"*" decay. A track used as a soft pion candidate must have a transverse momentum 
greater than 50 MeV, at least 25 wire hits measured by the tracking chamber, and a distance 
of closest approach to the event vertex smaller than 3 mm in the transverse plane. 

A cut on the transverse momentum of the D*^ vTg" system, > 1 GeV, is imposed in order 
to exclude the region of small acceptance of D*"*". 

The mass difference distribution AM = M(D°7r^) - M(D°) for the selected D° tt^ combi- 
nations in the two channels is shown in Figure ^ The contributions from the decay cascades 

(1) and (2) accumulate in narrow peaks close to the kinematic threshold. The mass difference 
spectrum is fitted by a sum of a Gaussian function for the signal and a term for the background 
of the form I]i=iaj(AM — m^)''% where Oj and bi are free parameters. The peak positions, 
determined by the fit, are 145.5 ± 0.2 MeV for the channel (1) and 146.1 ± 0.3 MeV for the 
channel (2), and agree well with the world average value for the mass difference mj)*+ —uij^o jllj . 
The good description of the background by the fit is corroborated by a background estimate 
obtained from the data themselves employing an event-mixing technique. For this, D*^ can- 
didates from a given event are combined with soft tracks from another event, containing D° 
candidates. The resulting background distributions are normalized to the data distributions in 
the region AM > 0.152 GeV and shown in Figure ^ by the dashed histograms. The number 
of reconstructed D*^ mesons is taken to be the number of observed entries in the signal region 
0.141 GeV < AM < 0.150 GeV, less the integral of the background fit component over that 
region. The 0*"*= signal is estimated to be 102 ± 17 events in channel (1) and 42 ± 11 in channel 

(2) . 

The combinatorial multiplicity in the signal regions AM < 0.150 GeV is 1.04 ± 0.01 for the 
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channel (1) and 1.05 ± 0.02 for the channel (2). There is no overlap of events in this region 
between the two channels and since the corresponding peak positions of the D*^ signal agree 
well, we add the two distributions shown in Figure C] and the resulting mass difference spectrum 
is shown in Figure |21 The total number of the observed D*^ mesons, obtained from the fit to 
the combined spectrum, is 149 ± 20. If the combined spectrum is split into two distributions 
for negative and positive charmed events, the fit result is 66 ± 14 D*~ mesons and 76 ± 15 D*"*" 
mesons. 



4 Inclusive D*^ Production Cross Section 

The cross section of inclusive D*^ production (summed over D*~ and D*"*") is determined in the 
visible kinematic region of experimental acceptance, to avoid model-dependent extrapolation 
uncertainties. In the present analysis, the selection cuts and the available statistics allow to 
cover the following phase space domain of D*^ pseudorapidity {rj^'l and transverse momentum 

D* 

Pt ■ 



\7]^'\ < 1.4 (3) 
1 GeV < p^' < 5 GeV. 

The differential spectra are obtained by fits to the mass difference distributions subdivided 
into three intervals of p^* or \r]^*\, the other variable being integrated over its kinematic region. 
Based on Monte Carlo studies, the resolution of the reconstructed p^* is determined to be about 
30 MeV and the resolution of \ri^* \ about 0.008 units of pseudorapidity. Thus the smearing and 
the resulting event migration between adjacent bins in the spectra of the reconstructed D*^ 
mesons is negligible. 

The efficiencies for the reconstruction of D*^ mesons are calculated separately for direct- 
photon processes and for single resolved-photon processes with Monte Carlo events generated by 
the PYTHIA program. A massive matrix element calculation with charmed quark mass value 
rric = 1.35 GeV and the SaSld parametrization of the parton distributions of the photon ^\ 
was used for the generation of events. The reconstruction efficiencies are calculated as a ratio 
of the combined number of reconstructed D*^ mesons in the two decay channels to the number 
of generated D*^ mesons and are presented in FigureElas functions of p^* and Evaluated 
in this way, the efficiencies take into account the corresponding branching fractions of the decay 
modes (1) and (2). The two sets of efficiencies are similar and agree within the errors. This 
implies that the relative proportion of direct and resolved-photon contributions to the charm 
production is not a major source of uncertainty in the determination of the D*^ differential 
cross sections in the phase space region defined by Equal contributions of both types of 
charm production processes in the kinematic region Q are assumed for the calculation of the 
reconstruction efficiencies, used for the cross-section evaluation. 

The measured cross sections of inclusive D*^ production, calculated as functions of p^* and 
|?7°*| and integrated over the corresponding bin, are listed in Table 1. When evaluating the 
differential cross sections, a correction based on the prediction of the combined Monte Carlo 
sample, used for the efficiency determination, is calculated such as to assign the differential 
cross sections to the centres of the corresponding bins. The values of this correction are found 
to be different from one only for the cross section set da/dp^ in the region of high p^* . The 
differential cross sections as a function of p^* , obtained after applying the bin-centre correction. 
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are listed in Table 1. The differential cross sections assigned to the bin centres are plotted in 
Figure Eland Figure El 



p^' \ GeVl 


A(Trr>pne \ pb 1 


da/dp^* [ pb / GeV] 


1/7°* 1 

1 / 1 


/S.a mf.no \ pb 1 


1-2 


92.9 ±22.2 ±16.7 


92.9 ±22.2 ±16.7 


0.0 - 0.4 


34.1 ±8.4 ±5.3 


2-3 


30.1 ±8.4 ±6.1 


28.0 ±7.8 ±5.7 


0.4-0.8 


47.5 ±11.0 ±9.6 


3-5 


11.3 ±3.9 ±3.0 


4.9 ± 1.7 ± 1.3 


0.8 - 1.4 


40.8 ± 15.8 ± 12.2 



Table 1: Measured cross sections Aameas for inclusive D*^ production, integrated over the 
corresponding bin. The third column of the table gives the differential cross sections as a 
function oip^* after bin-centre correction. The first errors are statistical, the second systematic. 



The systematic uncertainties on the measured cross sections are estimated by varying the 
selection of tracks and photons and by varying the cuts throughout the D*^ reconstruction. 
The contribution of the selection procedure to the systematic errors is in the range 8% — 17% 
affecting mostly the low-p^* region. The uncertainties in the K"??"*"?!" channel are higher than 
in the K^7r+ channel. The systematic uncertainties related to the background estimation are 
determined by using different forms for the background function in the mass difference fit 
and by changing the mass range of the fit and are found to vary from 5% to 10%. The D*^ 
reconstruction efficiencies are calculated also using a Monte Carlo sample generated by the 
JAMVG program which involves only direct processes in charm production, as well as for 
PYTHIA generations of direct and resolved processes with varied charmed quark mass value. 
The observed variations of the reconstruction efficiencies are taken into account as well as the 
Monte Carlo statistics, resulting in systematic changes of 5% — 14%. The contributions of the 
various sources of systematic errors are added in quadrature. 

The integrated cross section measured in the visible kinematic region is found to be 

a{e+e- e+e-D*±X ; 1 GeV < p?* < 5 GeV, \ri^* \ < 1.4) = 132 ± 22 ± 26 pb, 

where the first error is statistical and the second systematic. 

The integrated cross sections calculated separately for the K~7t~^ and K~7r+7r° channels, 
a = 124 ± 24 pb and cr = 142 ± 46 pb respectively (the errors are statistical only), agree well. 
This justifies combining the signals observed in the two decay modes. 

In Figure |3 and Figure El the differential cross sections are compared to next-to-leading 
order perturbative QCD computations, based on a massless approach in calculating the parton- 
level cross sections jTHj. In this scheme the charmed quark is considered to be one of the 
active flavours inside the photon. Three different sets of parton density parameterizations 
of the photon have been used in the calculations: GS [19^, AFG [20j and GRV-HO [21j. The 
renormalization scale, /xr, and the factorization scale of the photon structure function, /ip, have 
been taken as /xr = /xp/S = \l Pr + ml with charmed quark mass value mc = 1.5 GeV jTSj . 
There is a reasonable agreement between the data and the calculations. With regard to the 
variations of the predictions in the region of low transverse momentum, we should notice the 
limited applicability of the massless approach for px ~ rric f2'2\. 
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5 Summary 



The inclusive production of D* mesons in two-photon interactions at LEP is measured by 
reconstructing D*+ cascade decays involving D° decays into K~7t~^ and K~7r"'"7r° final states, 
as well as the charge conjugate decay chains. The integrated and differential cross sections 
of inclusive D*^ production are determined in the kinematic region 1 GeV < p^* < 5 GeV, 
1 77° I < 1.4 for which the acceptance is found to be insensitive to the relative mixture of direct 
and single resolved-photon processes. In this phase space domain the integrated cross section 
is measured to be a{e'^e~ — >• e+e^D*''=X) = 132 ± 22{stat.) ± 26{syst.) pb. A reasonable 
agreement is observed between the measured differential cross sections and the predictions 
based on next-to-leading order perturbative QCD calculations. 
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Figure 1: Mass difference distribution for D" decays into (a) K~7r~^ and (b) K~7r"'"7r°. The 
points are data, the hne is the resuh of the fit to the data points used to evaluate the D*"*" 
signal and the dashed histogram represents a background check, see the text. 
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Figure 2: Combined mass difference distribution for D'' decay channels K 7r+ and K 7r"'"7r°. 
The points are data and the hne is the result of the fit used to evaluate the D*+ signal. 
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Figure 3: Reconstruction efficiency of D*"^ mesons (including the branching fractions), deter- 
mined from PYTHIA generation of direct and single resolved-photon processes (a) as a function 
of p^* for < 1.4 and (b) as a function of |?7^*| for 1 GeV < p^* < 5 GeV. 
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Figure 4: The differential cross section of D*^ production as a function of tlie transverse 
momentum of tlie D*^ mesons for |?7^*| < 1.4. The points represent the data, the error bars 
show the statistical and systematic errors added in quadrature. The curves represent next-to- 
leading order QCD calculations ^H] for different parameterizations of the parton densities of 
the photon (GS [IH], AFG j^O] and GRV-HO [21j). 
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Figure 5: The differential cross section of D* production as a function of tlie pseudorapidity of 
the D*^ mesons for 1 GeV < < 5 GeV. The points represent the data, the error bars show 
the statistical and systematic errors added in quadrature. The curves represent next-to-leading 
order QCD calculations as in Figure El 
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